Entanglement of distant optomechanical systems 
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We theoretically investigate the possibility to generate non-classical states of optical and mechan- 
ical modes of optical cavities, distant from each other. A setup comprised of two identical cavities, 
each with one fixed and one movable mirror and coupled by an optical fiber, is studied in detail. We 
show that with such a setup there is potential to generate entanglement between the distant cavities, 
involving both optical and mechanical modes. The scheme is robust with respect to dissipation, and 
nonlocal correlations are found to exist in the steady state at finite temperatures. 
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I. INTRODUCTION 

Quantum entanglement is one of the most intriguing 
features of quantum mechanics. Although quantum me- 
chanics has proven to be highly successful in explaining 
physics at microscopic and subatomic scales, its validity 
at macroscopic or even mesoscopic scales is still debated. 
Some of the astonishing features appear when we try to 
apply quantum mechanical principles to macroscopic sys- 
tems. Superpositions of macroscopic systems is one ex- 
ample 

It is not yet completely clear to what extent quantum 
mechanics applies to macroscopic objects. Quantum phe- 
nomena such as entanglement generally do not appear in 
the macroscopic world. The difficulty of seeing quantum 
superpositions of macroscopic systems is often attributed 
to environment-induced decoherence. Such decoherence 
is thought to be the main cause reducing any quantum 
superposition to a classical statistical mixture [2j. Thus, 
an obvious but impractical choice would be to minimize 
the detrimental effect of the environment through per- 
fect isolation of the system of interest. Nonetheless, with 
the spectacular level of experimental advancements, the 
possibility of seeing macroscopic quantum superpositions 
appears to be within reach [3[. 

Related to this, quantum engineering Q in the field 
of optomechanics has made rapid advancement [j| . In a 
typical optomechanical setup, a mechanical system can 
be manipulated by radiation forces. Such systems have 
recently attracted much theoretical and experimental at- 
tention [f| 0]. This is partly because of their poten- 
tial usefulness in extremely sensitive sensor technology 
and in quantum information processing Q. Also, they 
are potentially one of the best tools to test fundamen- 
tals of quantum mechanics. Seminal progress has been 
made both theoretically and experimentally in this novel 
emerging field 0,0. 
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In a typical setting using optomechanical interaction, 
the main component is a cavity with a movable mirror. 
Light in the cavity and the movable mirror interact due to 
a coupling induced by the radiation pressure. As a result, 
the movable mirror executes harmonic motion around its 
equilibrium value 0, which thus alters the cavity reso- 
nance frequency. This, in turn, changes the circulating 
power in the cavity and hence the radiation pressure force 
acting on the movable mirror, leading to intrinsic nonlin- 
earities Strong light-matter coupling, both for opto- 
and for electro-mechanical systems, is a main ingredient 
in this emerging research field [1(J, LLl| ■ Within the strong 
coupling regime, radiation-pressure interaction has been 
successfully utilized for ground state cooling of mechan- 
ical oscillators (l^ - [T3 |. Some of the fascinating schemes 
include preparing the cavity mode and the movable mir- 
ror in a non-classical state [HI, [lH, preparing optome- 
chanical or fully mechanical Schrodinger cat states [l7l — 
[f9j . and even inducing quantum correlations between 
the subsystems j2fj| - |23j . Apart from the mostly studied 
cavity-movable mirror geometry, there have been some 
recent breakthroughs in exploring quantum features of 
a membrane in a cavity (24 . |25| . There are also recent 
proposals exploring the possibility of observing photonic 
analogs of the Josephson effect in an optomechanical set- 
ting (2<|. 

A common feature of most of these studies involve en- 
hancement of the radiation pressure coupling through in- 
tense laser driving of the cavity field. This is required 
to achieve strong radiation pressure coupling which oth- 
erwise is too weak to observe any non-classical phe- 
nomenon. Although most of these studies are restricted 
to Gaussian state preparation involving optomechani- 
cal interaction, there have been some recent proposals 
to study non-Gaussian quantum states in the regime of 
single-photon optomechanics p7l . [28j . 

Motivated by these theoretical and experimental ad- 
vancements, we shall here explore the possibility of entan- 
gling mechanical and optical modes of two distant cavi- 
ties. In previous studies of the entanglement of distant 
mechanical mirrors, squeezed light was used as an avail- 
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able resource in order to entangle the two distant mirrors, 
which were either part of the same optical cavity [HI, [2l[ 
or belonged to two different cavities [l^]. In the present 
work we are interested in a physically different setup, 
where two distant Fabry-Perot cavities each fabricated 
with one movable mirror are coupled by an optical fiber. 
We show that, as a result of a combination of the op- 
tomechanical interaction and an optical-fiber mediated 
coupling, the two distant optical and mechanical modes 
become entangled. Moreover, we explicitly study two 
different regimes of physical interest. First we will im- 
pose an approximation of Born-Oppenheimer (BO) type, 
and consider a scenario in which the two cavities are not 
externally pumped. In this regime the two mechanical 
modes are found not to be very strongly entangled. The 
advantage is however that the approximation allows us to 
find an analytical solution describing the evolution of the 
state of the two mirrors. Thereafter we work in a regime 
where the coupled cavities are strongly driven. Here we 
explicitly derive the relevant quantum Langevin equa- 
tions (QLE) and construct the covariance matrix gov- 
erning the dynamics of all the optical and mechanical 
modes. 

The paper is organized in the following way. Section [TT1 
introduces the theoretical model and the physical setup. 
This forms the framework for Sec. MI Bl in which the 
unitary evolution of the system is studied, followed by 
analysis of the dissipative regime in Sec. IIII CI and a brief 
discussion on the validity of the BO approximation in 
Sec. MIDI A quantum Langevin approach is introduced 
in Sec. IIV1 and we conclude with a short discussion in 
Sec. El 



II. PHYSICAL SETUP 

We consider a physical setup comprised of two identical 
Fabry-Perot cavities, each with one fixed and one mov- 
able mirror, as schematically shown in Fig.[T] We assume 
that only one resonant mode of each cavity is populated, 
and that these two modes arc coupled via an optical fiber. 
The two modes have the same frequency, lo = 2wc/L; 
where L is the cavity length, and are described by the 
creation (annihilation) operators a^(a) and b'(b), respec- 
tively. Furthermore, we assume that each movable mir- 
ror has been cooled near to its ground state, so that it 
is operating in the quantum regime. Under the action 
of cavity-photon-induced radiation pressure, the movable 
mirrors will oscillate about their equilibrium positions. 

If we assume that the two mirrors move distances x and 
y along the respective cavity axes, so that the two dis- 
placements are much smaller than the wavelength of each 
cavity mode in one cavity round-trip time, then scatter- 
ing of photons to other cavity modes can be safely ne- 
glected [H, [30| . The effective lengths of the cavities will 
then become L + x and L + y, with new resonance fre- 
quencies u a = 2irc/(L + x) and u>b = 2nc/ (L + y), where 
x and y are the instantaneous displacements of the two 
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FIG. 1: (Color online) Sketch of the physical setup to en- 
tangle distant optomechanical modes. Two optomechanical 
cavities pumped by classical laser fields are coupled to each 
other by an optical fibre. As a result of indirect coupling 
mediated by the two cavity modes, the two movable mirrors 
become entangled. Furthermore, two initially uncorrelated 
auxiliary cavity modes interact independently with the two 
entangled movable mirrors, which induces non-local corre- 
lations between the two modes. Using standard homodyne 
measurement techniques non-local correlations between the 
two auxiliary cavity modes can be read out giving an indirect 
signature of quantum correlations between the two mirrors. 



cavity mirrors from their equilibrium positions. With the 
above assumption, i.e. x/L,y/L <C 1, the free evolution 
of the two optical cavity modes in the adiabatic regime 
takes the form 1301 
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Under the action of a weak radiation-pressure force, each 
movable mirror undergoes small-amplitude oscillations 
with frequency tt. In the absence of external driving, 
the full Hamiltonian of the two coupled cavities thus be- 
comes 



H = Hf, 



^ y2 + i^ + 
2m 2 2m 



(2) 

where A is the inter-mode coupling between the two cavi- 
ties. This coupling could be mediated by, e.g., an optical 
fiber connecting the two distant cavities. Introducing 
dimensionless conjugate variables qi and pi for the ith 
movable mirror, ([2]) can be rewritten as 



H 



Ml 



hw {a} a + tfbj + ^ (gi 2 + p\ 2 ) + ^ {q 2 + V~i) 
+h\ (afb + tfa^J - hg (a^aq'x + Pbq-^j 



(3) 



where g = (lu / ' L)^/h/m^l is the radiation pressure- 
induced coupling between the cavity modes and the mov- 
able mirrors. The Hamiltonian ([3]) will form the basis for 
the analysis in the next section, where we will make an 
adiabatic approximation that allows us to study the uni- 
tary evolution of the two movable mirrors. 
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III. BORN-OPPENHEIMER APPROACH 
A. Effective adiabatic model 

The frequency mismatch between optical (uj/2tt ~ 
10 14 Hz) and mechanical (ft/27r - 10 6 - 10 9 Hz) de- 
grees of freedom is enormous [3] . This suggests a separa- 
tion of the Hamiltonian ([3]) into two parts, one with very 
rapidly evolving optical modes and another with slowly 
varying mechanical modes. In the limit that the mirror 
coordinates qi and q\ remain stationary with respect to 
the rapidly evolving cavity modes a and b, we can diag- 
onalizc the interaction between the two cavity modes of 
Hamiltonian ([3]). 

We first introduce the collective excitation operators 
A and B obeying 



two mirrors evolve, 



co&6 smB 
— sin# cos# 



(4) 



Choosing tan 29 = 2\/{g(q\ — q 2 )) and substituting for 
the new field modes, the rapidly varying optical part of 
Hamiltonian (J3|) reduces to 



h uj 



qi + <?2 



The Hamiltonian describing the dynamics of the two 
movable cavity mirrors thus takes the form 



flmir = fi— (q\ 2 



~ 2 , » 2 , - 2> 
Pi +Q2 +P2 j 



(0) 



The treatment this far is exact. Typically, the cavity 
field will adiabatically follow the slow motion of the two 
movable mirrors. Thus by considering rapidly varying 
collective cavity modes and slowly var ying m irror modes 
we can make the BO approximation [3lL l32j , and write 
the collective wave function of the cavity-mirror coupled 
system at time t as 



|¥(t)> = J2 P(n A ,n B )\n)\$(n,t)). 



(7) 



Here, P{n A ,n B ) is the probability distribution of the 
collective cavity fields, \n) = \n A ,n B ) denotes the time- 
independent index of the energy levels of the two col- 
lective cavity modes A and B, in the adiabatic limit, 
in which A^A|n) = n A \n ) and B^Bl n) = n B \n), and 

= e _iHmir */ ri |$(n, 0)) denotes the time-evolved 
wave function of the two movable mirrors [33|, HH . The 
approximation in assigning a system wave function of the 
form ([7]) lies in the fact that the coefficients P(n A ,n B ) 
are time-independent, and as a consequence no popula- 
tion transfer occurs between different photon states |n). 

Within this BO approximation, the cavity modes can 
be seen as inducing an effective potential in which the 



- .9 qi 77 92 ) A + n B ) 
-h^g^ qi - q2 y+AV nA ~ nB . (8) 



Since we have assumed that the oscillation amplitudes 
of the movable mirrors are small, it follows that their 
relative displacement q\ — q 2 is also small. Therefore, 
it is sufficient to expand the second term in the cavity 
Hamiltonian (JSJ) to second (quadratic) order in q\ — q 2 . 
This can be justified, since for a typical optomechanical 
cavity with optical frequency uj/2tt ~ 10 14 Hz, length 
L ~ 1 mm, mirror frequency Q,/2tt ~ 10 6 Hz, and 
with a zero-point-oscillation amplitude of 0.02 pm, one 
finds that g ~ I0 4 Hz. With the reasonable estimate 
A = 10 5 Hz one gets (g/2\) 2 ~ 10" 3 . This results in 
an effective Hamiltonian describing the dynamics of two 
coupled movable mirror in absence of any losses, 



H, 



, fc / - 2 , - 2 , -2, -2-1 
r ~ ' Ql Pl q2 P2 > 



-(n A - n B)hX— 2 (qi - q 2 ) 2 , 



(9) 



where we have dropped all the constant and classical 
driving terms from the Hamiltonian. Dynamical prop- 
erties of entanglement in a model related to the one of 
Eq. © was recently studied for a closed system f35j . 
Equation ((6]) can be rewritten in terms of 
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such that 

H mil = HQ (c^c + Sctj (10) 
-(n A - n B )h\ (j^J (c 2 + £t 2 + 2#c) 

-(n A - n B )h\ (j^j (d 2 + d) 2 + 2d)d) 



+{n A -n B )h\ ( ^ ) (c + c f ) (d + S 
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Introducing center-of-mass and relative modes, 



Eq. (|T0|) becomes 

H mir = hflC^C + h(Q- 4N\)DW (12) 

-2ivriA (p 2 + b ]2 
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where TV = (ha — ns)(<7/4A) 2 . The Hamiltonian in the 
above form is known to generate squeezing in the D- 
mode [3a |. which will also be manifested as quantum 
correlations between the two mirror oscillations. 

After arriving at this simplified form of the Hamilto- 
nian governing the dynamics of the two movable mirrors, 
we now provide a fully analytical treatment describing 
the state evolution of the two mirrors. Firstly we shall 
discuss the unitary dynamics of the system in section 
IIIIBI and provide a closed-form expression for the time- 
evolved mirror operators c(t) and d(t) in the Heisenberg 
picture. This will allow us to solve for the dynamics of 
initially uncoupled movable mirrors for an arbitrary ini- 
tial state. This will then be followed by IIII CI where we 
shall provide a full solution of the master equation de- 
scribing the dissipative dynamics of the two indirectly 
coupled movable mirrors. 



giving the closed expressions for the time evolved opera- 
tors c(t) and d(t), 

c(t) = - F(t)c(0) + G(t)d{0)+2isin(2u> t)uvc\0) 

-2iam(2u) Q t)uvd f (0) 
d(t) = - G(t)c(0) + F(t)d{0)-2isin(2u> t)uvc\0) 



+2ism(2ui t)uvd^(0) 



(18) 



where F(t) and G(t) are time-dependent complex func- 
tions given by 



F(t): 

G(t) 



e -iUt + u 2 e -i2«, t 
. e -iQt + y 2 e i2uj t _ 



— v 2 e i2LJot 



(19) 

(20) 



B. Unitary evolution 

The Hamiltonian (|T2"j) describing the dynamics of the 
two movable mirrors can be further diagonalized by a 
Bogoliubov transformation. We define E and E such 
that 



D J ~ \ v u ) \ E 



(13) 



where u 2 — v 2 — 1. Setting 




(14) 

the Hamiltonian (|T2"j) reduces to the diagonal form 

H mir = htt&C + 2fiw &E, (15) 



where 



LU 



A I 



(n-8XN)fl 



2^1~AM 2 (n-A\N)' 
2NX 



Q-4NX 



(16) 



We can then straightforwardly solve the equations of mo- 
tion for the operators C(t) and E(t), 



C(t) = C(0)e iM 
E(t) = E(0)e- l2uJot , 



(17) 



With the solution of the operators c(t) and d(t) now 
in hand we can faithfully describe the unitary dynamics 
of the two movable mirrors for any arbitrary initial state. 
Of particular interest are initial Gaussian states includ- 
ing thermal, coherent, and squeezed states. A Gaussian 
continuous variable state can be fully described in terms 
of a real symmetric covariance matrix V. For a two-mode 
Gaussian continuous variable system, the covariance ma- 
trix V can be written as 



V = 



A C 
C T B 



(21) 



where T denotes matrix transpose, 
A 



{(c + c)) 2 )/2 ([c + ct,i(ct-c)] + )/2 
* .Ytf_g)] + )/2 <(i(ct -c)) 2 )/2 



([c + c\i{c 



B = 



((«2+#) a }/2 ([d + - d)}+)/2 

([d + rff,i(rff - d)] + )/2 - d)) 2 )/2 

((c + #)(d + S))/2 (i(c+ct)(dt-d))/2 \ 
■ ' k ((i(gt - c)(d + S))/2 (-(ct - c)(dt - d))/2 ) ' 



and ([fi,fj} + ) = {{fifj +fjfi)))/2. Once we have the 
covariance matrix, all the quantum statistical proper- 
ties of Gaussian continuous variable states can be con- 
structed. Also worth mentioning is the important fact 
that the Hamiltonian (|10[) is quadratic in the position 
and momentum coordinates of the movable mirrors. An 
initial Gaussian state of the mirror evolving under (|10[) 
will therefore maintain its Gaussian character. 

A widely used entanglement measure is the logarith- 
mic negativity, which is an entanglement monotone and 
fairly easy to compute [37j |. For a two- mode Gaussian 
continuous variable state characterized by the covariance 
matrix V, logarithmic negativity is defined as 



J V = Max[0,log(2f>_)], 



(22) 
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FIG. 2: (Color online) Time evolution of the degree of en- 
tanglement, as measured by the logarithmic negativity, as a 
function of initial temperature of the movable mirrors, mea- 
sured in terms of n t hermai- The dimensionless parameters are 
chosen such that O = 1, g = 10~ 2 , A = 10 _1 , otA = 4 and 
olb = 1. 



where V- is the smallest of the symplectic eigenvalues of 
the partially transposed covariance matrix given by 



= J a/2 - ^J{o 2 -4DctV)/2 



cr = DetA + DctB - 2DetC. 



(23) 



We analytically reconstruct the time-dependent covari- 
ance matrix, from which it is then straightforward to 
compute the logarithmic negativity, with a typical so- 
lution shown in Fig. [2j In these calculations, the log- 
arithmic negativity has been weighted with a coherent 
state probability distribution for the collective cavity 
modes A and B such that P(n A ,n B ) = exp(— (la^l 2 + 
\ot B \ 2 ))\a A \ 2nA \a B \ 2nB /(n A \n B \). Such averaging ac- 
counts for initial quantum fluctuations in the two cavity 
modes. A non-zero value of M quantifies the degree of en- 
tanglement between the two movable mirrors. As can be 
seen from Fig. [2J increasing the initial temperature of the 
mirror degrades the quantum correlations and eventually 
leads to completely separable states of the two mirrors. 
The figure also gives a clear example of entanglement 
sudden death and birth [38[, arising from the common 
coupling of the mirrors to the two cavity modes. 



C. Dissipative dynamics 

In any physical setting, coupling to the environment 
is inevitable and typically results in decoherence of the 
quantum state to its classical counterpart. In the scheme 
of interest to us, there can be two main causes of dissipa- 
tion. One is the photon leakage through the two cavities, 
and the other is thermal decay of the states of the mov- 



able mirrors due to their coupling to baths of non-zero 
temperature. 

So far we have treated the cavity modes in the adia- 
batic approximation by expanding the global wave func- 
tion in terms of energy eigenstates of the collective cavity 
operators A and B. Although it might look somewhat 
artificial at first, this approach has its own advantages. 
Firstly, it allows us to derive a closed-form analytical re- 
sult governing the dynamics of the two coupled mirrors, 
which is otherwise a nonlinear problem in itself. Sec- 
ondly, it also provides us with useful physical insights into 
how the cavity assisted entanglement of spatially sepa- 
rated mechanical oscillators originates. In the present 
section we will continue treating the two cavity modes 
in this semiclassical regime and defer the explicit calcu- 
lations involving reservoir induced quantum fluctuations 
in the cavity modes till the next section. More precisely, 
we assign coherent states for the two cavity modes and in- 
troduce cavity losses in terms of a non-Hcrmitian Hamil- 
tonian. 

A phcnomcnological way to introduce cavity losses is 
to shift the cavity resonance frequency u by —in where k 
is the cavity decay rate. Then, under the BO approxima- 
tion, the two indirectly coupled movable mirrors evolve 
according to the non-Hermitian Hamiltonian 



in(n A + n B ) 



(24) 



where H mir is given by (fTU)) and we again have neglected 
all the classical driving terms in the Hamiltonian. Apart 
from the cavity losses, the two cavity mirrors might un- 
dergo further decoherence due to their inevitable cou- 
pling to the external environment. The time evolution of 
the mixed state of the two movable mirrors obtained by 
tracing over the cavity field distribution takes the form 



Pmir(i) 



I 



E 



,(*) 



(25) 



where 

IV n, a .n. 



,(*) 



cxp I 



\a A e- Kt \ 2 ) \a A e- Kt \ 2nA 



xexp(-|ase K '| 2 ) \a B e 



(26) 
%B /n A \n B \, 



and /3 lnir (i) is the time-evolved reduced density matrix of 
the two movable mirrors with the photon number differ- 
ence n = riA — n B . It turns out that if both collective 
cavity modes are initially coherent states, this particu- 
lar method of taking dissipation into account is not only 
accurate but also exact [39j . This is because an initial 
coherent state evolving in a purely dissipative channel re- 
mains a coherent state, although with an exponentially 
decaying amplitude [40| . The time evolution of pj^ T (t) 
in the Born-Markov approximation is then described by 
the Lindblad-type master equation [4l|, [42[ 



mil") pmix 



-nC c \p 



r ~r ^ Tl^d^ Pmir 



r r 

+ -(n + l)£ ci o rair + — (n + l)£ d /5mir, (27) 
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where r is the decay rate of each movable mirror due to 
its coupling to a heat bath with average thermal occu- 
pancy n, and C x p = 2xpx' — x'xp — px^x. In terms of 
the center-of-mass mode C and relative mode D, Eq. ([27)1 
can be equivalently written as 



dp„ 



dt 



r r 

+ -(n + l)£c/5mir + g (" + iKDPmir, (28) 



where -ff m i r is given by fj 12|) . 

To solve the master equation (|28|). we define the 
normal-ordered quantum characteristic function [ill . 
|42T | for the two movable mirrors as x(e, r],t) = 
(e^e -6 CgyD^ e -y £>y Using standard quantum opti- 
cal techniques [ill, |42| , the master equation (|28| can be 
rewritten as a partial differential equation for the quan- 
tum characteristic function x(e, rj, t) of the form 



d_ 

di 
where 



X (e,V,t) ^ z T M\/ X (e,r),t) + A\z T Kz X (e,v,t), (29) 
z T = (u 1 ,u 2 ,v 1 ,v 2 ), 



V = ( 9 

\ dui ' du2 5 dvi ' 9i)2 



"2 



''1 



t'2 



<ic + £d + (- c +<i d 

i2s/2 

ec-gd+e*-gj 
2V2 

e c -e d -e;+e; 
i2V2 



(30) 



and 



e + f] 



The 4x4 matrix coefficients of Eq. ([29]) read 



V 



M 



with 
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FIG. 3: (Color online) Temporal evolution of the degree of 
entanglement between two indirectly coupled movable mirrors 
as measured by the logarithmic negativity. Compared with 
Fig. [2l here losses in all modes have been considered and 
the degree of entanglement is consequently somewhat smaller, 
but importantly, it survives for a reasonably long time. Each 
mirror is initially assumed to be in its ground state and the 
dimensionless parameters used are chosen such that Q — 1, 



9 = 10-, 

and n — 



A = 10 , (xa = 4, a B = 1, k = 10" 



T = 10" 



For an initial Gaussian state of the two movable mir- 
rors, it is consistent to make the following Ansatz for the 
quantum characteristic function, 



X(e, f), t) = exp [— z T L(t) z + iz T q(t)] 



(34) 



where L(t) is a 4 x 4 time-dependent symmetric matrix 
and q(t) is a 4 x 1 time-dependent vector. Using the 
above Ansatz in Eq. (f2U)) results in the coupled matrix 
differential equations 



L = ML + LM 

q = Mq. 



4AK, (35) 
(36) 

Making use of the fact that L is a 4 x 4 symmetric matrix, 
it can be decomposed into 2x2 square matrices such that 



L -( P Q 
VQ T R 



(37) 



where P and R are 2x2 symmetric matrices. Obtaining 
an explicit form for the time-dependent quantum char- 
acteristic function \(e,i],t) now reduces to solving 2x2 
coupled matrix differential equations. 

Although an exact analytical solution can be arrived 
at, it is too lengthy to be reported here. Nonetheless, 
the time-dependent covariancc matrix V can be fully 
reconstructed from the quantum characteristic function 
x(e,7y,i). This can be easily seen by noting that from 
the quantum characteristic function one can obtain the 
expectation values of quantum mechanical observables, 
e.g., (&™(t)$»(t)) = (^-)-(^-)» x ( ec , e(2 , *)l fc ,e d =0 
and thus all the elements of the covariance matrix can 
be found. 

As a measure of entanglement between the distant cav- 
ity mirrors we again compute the logarithmic negativity. 
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The result of such a calculation is shown in Fig. [3] As is 
clear from the figure, under the action of cavity-mediated 
coupling, the two movable mirrors exhibit entanglement. 
Although the entanglement generated is not too large, 
it is sustained over a reasonably long timescale. The de- 
gree of inseparability between the two mirrors can be im- 
proved significantly either by a conditional measurement 
of the cavity field or by increasing the difference in the 
mean number of photons in the field distributions of the 
two cavity modes. Thus we conclude that the aforemen- 
tioned protocol is indeed capable of generating quantum 
entangled states of two movable mirrors, which are robust 
with respect to dissipation for a long time. It should be 
pointed out that the logarithmic negativity approaches 
zero exponentially for large times due to the decay of 
photons out of the cavities. In order to have sustainable 
non- vanishing entanglement, the photon modes must be 
driven externally to prevent the absence of photons. This 
will be discussed in the following section. 



looked so far concerns the fact that we consider an open 
quantum system for which the concept of adiabaticity 
must be handled with care [50]. Loosely speaking, one 
could imagine non-adiabatic transitions induced by the 
reservoir. In the following discussion we will argue that 
such transitions are indeed not present in our model. 

Recently an extension of the BO approximation has 
been presented for a quantum system coupled to a large 
reservoir j5fj| . Following this work, we can write the 
Hamiltonian of the closed quantum system as H = 
H S (X) + H f (X,Y), where H S (X) is the Hamiltonian of 
the slowly varying dynamical variable X and Hf(X,Y) 
is the interaction Hamiltonian between the slowly vary- 
ing variable and the quickly varying dynamical variable 
Y. If the coupling between the quantum system and the 
reservoir is such that its evolution can be described using 
the Lindblad approach, then, in the case of no quantum 
jumps, the evolution of the quantum system is governed 
by the equation 



D. Validity of the Born-Oppenheimer 
approximation 

So far we have used the BO approximation to separate 
the slow dynamics of the two movable mirrors from the 
rapidly evolving population of the two cavity modes. The 
BO approach has proved to be a fundamental tool in 
various applications of physics and quantum chemistry 
[3ll l32j . It is nothing but an extension of the quantum 
adiabatic theorem to a quantum system with two sets of 
variables whose dynamics can be separated due to very 
different dynamical time scales. 

The original version of the adiabatic theorem asserts 
that an initial eigenstate of a slowly varying Hamilto- 
nian will remain in the same instantaneous eigenstate 
throughout the evolution. In our system, the absence 
of transitions between instantaneous eigenstates trans- 
lates to fixed probabilities for the states |n), see Eq. (J7J, 
to be occupied. There has been continued interest in 
the application of the adiabatic theorem in slowly evolv- 
ing quantum mechanical systems including Berry's phase 
and effective gau ge theories p3l |44| . geometric quan- 
tum computation |45j , and adiabatic quantum computa- 
tion [46|. However, recently the generally accepted crite- 
ria for adiabaticity in quantum physics has suffered some 
criticisms [47|, which has lead to various arguments and 
counter-arguments to ascertain the justification of the 
adiabatic approximation in quantum theory (48l . |49| . A 
breakdown of adiabaticity seems to appear only in rather 
special cases of rapid resonant driving of the system [4!| . 
This kind of external resonant driving is absent in our 
model and we expect a standard application of the adia- 
batic theorem, or more precisely an application of the BO 
approximation, to be unproblematic. In the beginning of 
Sec. IIII Al wc pointed out that the two types of oscilla- 
tors evolve on very different time scales which normally 
guarantees adiabatic dynamics. However, one issue over- 



(38) 



Here H off = H S (X) + H' f (X,Y), and H' f {X,Y) is the 
non-hermitian part of the Hamiltonian, with dissipation 
assumed to affect only the quickly varying variable of the 
Hamiltonian. Treating the slowly varying variable X as 
a parameter, we can solve for left and right eigenstates 
(^nx(X)\ an d |*nxP0) 01 tnc non-hermitian Hamilto- 
nian H'f{X, Y), with complex eigenvalues E n> x(Y). Ex- 
panding the eigenstates of H e s in terms of xO^)) wc 
get 



N 



|$> = £c„K(X)>K >x (Y)>, 



(39) 



where N is the dimension of the quickly varying vari- 
able Y and c„ are the expansion coefficients. It can be 
shown that the eigenvalue equation H e g\&) = E\$>) can 
be recast in the form {H a + Hp)tp = Etp [5fJ, where 



Hx+ExiY) 

H 2 + E 2 (Y) 






H N + E N (Y) 



H (X) 



H P (X) = 



Here H n = H n>n , E n {Y) = E n , x {Y) and H n . m {X) = 
(*^(y)|i? s (l)|^ x (F)). Treating H P as a small per- 
turbation and using standard time independent pertur- 
bation theory, corrections to all orders for the eigenstates 
and the eigenenergies can be calculated. Since there is 
no resonant driving, we can conclude that the BO ap- 
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proximation is accurate as long as the condition 
<< 1 



<*»i^.»i*3?> 



E 



[o] 



E 



[0] 



for all fc , n ^= k,n. 

(40) 

is satisfied |5(j |. 

Returning to consider our original physical system of 
two coupled optomechanical cavities, the non-hermitian 
Hamiltonian containing the cavity decay terms takes the 
form 



He, 



= h [ ui — in — g 



qi + 92 



(A^A + B^B) 



-hy 9 \q1-q2) +4A 2 



+^(q 2 i + Pi + ql+pl)- 

The quickly varying part of the Hamiltonian 

H f = h[u J -iK-g{q 1 +q 2 )/2]{A^A + B^B) 



(41) 



hyg^qi - q-2? + 4A 2 (iU _ &B)/2 

has left and right eigenstates ($^x(^)l = { n A\{nB\ an d 
|f^ x (y)) = \nA)\nB) respectively, with complex eigen- 
values given by 



E n ,x(Y) = h(u} - in - g(qi + q 2 )/2)(n A + n B ) 



Wg 2 (qi -q?) 2 + 4A 2 (tm - n B )/2. 

As a passing remark, the perturbation terms Hp are in 
general non-zero. This demands that the dissipation in 
the slowly varying variables must be slower than the dissi- 
pation in the quickly varying variables. If this condition 
is not met, then strong dissipation in the slowly vary- 
ing variables enlarges the perturbation Hp, which might 
make the BO approximation invalid. Roughly speak- 
ing, a stronger rate of dissipation in the slowly varying 
variables accelerates their rate of change, thus making 
it hard to distinguish their dynamics from that of the 
quickly varying variables. The slowly varying variables 
then evolve under the effective Hamiltonian 



^ (q 2 + pi + q 2 + pl) 



W5 2 (9i-92)' + 4A 2 



(n A - n B ) 



-h ( u — ik — g ^ ) {riA + tib)- (42) 



It should be noticed that in this treatment, the off- 
diagonal terms of the perturbation Hp , which are given 
by H rhm {X) = ( 

TlA I (j^B I -^mirror |toa)|?ti_b), are identi- 
cally zero. This implies that the zeroth order BO ap- 
proximation for this particular physical model is not only 
accurate but also exact. The conclusion remains the same 
even when dissipation in the slow variables is taken into 
account 1501] . 



IV. QUANTUM LANGEVIN APPROACH 



As mentioned before, in a typical optomechanical 
setup, due to coupling induced by radiation pressure, a 
movable mirror interacts with a cavity mode. Unfortu- 
nately, the radiation pressure coupling for an undriven 
cavity with a movable mirror is usually very weak. This 
problem can be circumvented by driving the cavity with 
a coherent classical laser field. Driving with an intense 
laser field enhances the radiation pressure coupling and 
thus facilitates the observation of non-classical phenom- 
ena such as entanglement between mechanical oscillators 
and light. 

In the previous section, cavity driving was not taken 
into account. Instead the two cavity fields were assumed 
to be initially in coherent states and the system time- 
evolution in the presence of the decay of the two cav- 
ity fields was studied. In what follows, we shall instead 
study the situation of cavity driving and show that ro- 
bust steady state entanglement may exist between differ- 
ent optical and mechanical modes. To this end, we find it 
more convenient to work in the Heisenberg picture. For 
our system of two coupled cavities with movable mir- 
rors, the Hamiltonian in the interaction picture of the 
two driving lasers with frequency ujl now takes the form 



H 

7 



o 



(u - oj L ){a)a + Pb) + -{q\ +p\ 2 ) + + pl) 



+X(a ji b + tfa) — ga* aq\ — gb^bq2 
+irj{a i — a) + irj{W — b). 



(43) 



Here 77 = ^/2P c «;/fiw c is related to the driving laser, 
where P c is the power of the driving laser and n is the 
damping rate, identical for both cavities. 

The Hamiltonian (|43|) describes the closed-system dy- 
namics of the two driven coupled cavities with movable 
mirrors. However, as discussed in the previous section, 
the dynamics of the system is also affected by damp- 
ing and noise. The main channels of dissipation in our 
system are the decay in the cavity modes and the cou- 
pling of the movable mirrors to their independent ther- 
mal baths [5l| . One possible way to take into account 
all the damping and noise processes is to use quantum 
Langevin equations (QLEs). The QLEs are equivalent to 
the Heisenberg equations of motion for time-evolving op- 
erators, where noise and dissipative processes have been 
included phenomenologically |4l| . 

For the Hamiltonian (|4"3")l . the QLEs for the cavity and 
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the mirror modes become 



the optical and mechanical modes, we get 



da 
dt 

db 
dt ' 

dpi 
dt 

dp2 
dt 



dt 

dq-2 
dt 



—i\b + igqid + rj — (k + iA)a + 
—iXa + igq 2 b + 77 — (k + iA)b + 
-f2gi + go? a - -f m pi + ii{t), 
-VLq 2 + gtib - j m p 2 + i 2 (t), 

Qp2, 



Ik}) j 7r . 



(44) 



where A = u> — ujl is the laser detuning from the cavity 
resonance frequency w, n is the decay rate of each cav- 
ity and 7 m is the thermal decay rate, identical for the 
two movable mirrors subjected to independent Brownian 
motion noise characterized by the operators ix(t) and 
i 2 (t) respectively. The quantum noise operators have 
the quantum statistical properties 



Qi = 



g\g s \ n 



1 2 — ' 



Pf = 0, 



_ —i\ri+i){K+iA b ) 
s ~ A 2 +K 2 +iK(A a +A b )-A a A b ' 



bs = 



— i\ri+rj{K+i/S. a ) 



A 2 +K 2 +iK(A a +A(,)-A Q A [ , ' 

A a = A - gq{ = w - oj l - gq{, 
A5 = A - gq s 2 = uj - uj l - gq s 2 . 



(47) 



(h(t)) = (e 2 (t)) = 

mm 1 )) = ^ f e- iJ[t 



(45) 



1 + coth 



hjj 



5ij, 



where i, j £ 1, 2 and Sij is the Kronecker delta. We 
have also introduced independent cavity input noise op- 
erators, di n (t) and bi n (t) for the first and second cavity 
respectively. For the case of optical fields, huj/ksT >• 1, 
and hence the mean number of thermal photons can be 
safely neglected. In this limit the noise operators di n (t) 
and bi n (t) satisfy the two-time correlations 



(a in (t)al(i/))=6(t-t'), 
(al(t)a in (f))=0, 
(kn(t)bl(t')) = 8(t-t'), 
(bl(t)hn(t')) =0. 



(46) 



We arc interested in investigating the possibility of 
achieving steady-state entanglement between distant op- 
tical and mechanical modes. To pursue this aim we have 
to solve the set of coupled nonlinear QLEs (|4"4"]) . This task 
is difficult, but it is simplified in the presence of strong 
external driving, in which case linearization of the set of 
QLEs around the steady-state values is justified. Solving 
the set of QLEs (|4"4")l for the steady-state amplitudes of 



In the regime where the two cavities are very intensely 
driven, such that |a s |,|6 s | ^ 1, and by expanding the 
mechanical and optical mode operators as quantum fluc- 
tuations around the steady state values (a = a s + 5d, 
b = b s + 5b, fa = qf + Scji and f>i = pf + 5p l for i = 1, 2) 
we obtain the following linearized QLEs for the quantum 
fluctuations, 



dSd 
dt 

dSb 
dt 

dSqi 
~dT 



dt 

dSpi 
dt 



dt 



—8d(n + iA a ) — i\5b + iga s 5x + \/2nai n . 
—Sb(n + iA ) — i\Sa + igb s 8y + y/2nbi n , 

= rwp 2 , 

= —QSqi + g (\a s \ 2 + a* s Sa + a s 5a^) , 
= -tl5q 2 + g (\b s \ 2 + b* s 5"b + b s 5b^ 



(48) 



By further introducing the position and momentum 
quadratures for the two cavity modes and their input 
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noises, 



d6Xg 
dt 

dSPa 

dt 
dSX b 

dt 
dSP b 
dt 
dSX9 n 
dt 

dt 
dSXl 

dt 
d8X b in 

dt 



we can rewrite Eq. 



d(6tf + Sd) 
dt ' 
.d(Sa) — Sd) 
1 dt ' 
d(6# + 5b) 
dt ' 
d{5tf -5b) 
1 dt ' 

d (8dl n + Sajn) 
dt 

■ d (8d\ n - 5d in ) 

1 dt 

d(Sbl + 5b in ) 
dt 

d(Sbl-Sb in ) 
dt 

in the compact form 



d 4 = ZR + N. 

dt 



(49) 



(50) 



Here 

R T = (Sqi,Spi,Sq2 1 Sp2 1 SX ai SP ai SX bl SP b ), 
N T = (0,^(^,0, e 2 (t),V2HSX^ n (t), 

V^5P? n (t) , V^SXl (t) , V2HSPI (t)) , 

and 



Zi Z 2 
Z2 Z3 



(51) 



with the 4x4 matrices 



Zi = 














-0 


-7m 

















n 








-n 


~7m 



z, = 





g s a 



g s b 

-k A a A 

-A a — k -A 

A -k A b 

-A -A b -k 



(52) 



The phase reference has been chosen such that a s and b s 
are real with g s a = V2ga s and gl = s/2gb s . In Eq. (|5U|) . 
Z is the drift matrix. Stability (in the steady state) de- 
mands that the real part of all the eigenvalues of Z must 




FIG. 4: (Color online) Logarithmic negativity as a measure 
of entanglement between (a) two distant cavity mirrors, (b) a 
mirror and adjacent cavity mode, and (c) a mirror and distant 
cavity mode, plotted as a function of detuning A and average 
thermal occupancy of the two mirrors fii = ni = n. We have 
chosen the different physical parameters such that fl = 1, 
g s a = g s b = 2.5, A = 20, k = 0.08, 7m = 0.01, and A a = A b = 
A 



be negative. All the parameters in the present work have 
been chosen such that the system is stable in the steady 
state. 

The dynamics of the coupled cavities with movable 
mirrors is governed by the first order matrix differen- 
tial equation (foT)|) . For an initial Gaussian state of the 
two cavities and their movable mirrors it is sufficient to 
fully characterize all the quantum correlations by explic- 
itly evaluating the 8x8 symmetric covariance matrix 
V where Vij(t) = ((R t (t)Rj(t) + R j (t)R l (t)))/2. If the 
system is stable in the steady state, then the covariance 
matrix takes the form 



V, 



E 



ds I ds'W^Wj.gis'^p.gis-s 1 ), 
'0 



(53) 

where W = cxp(Zs) and $ p , q (s - s') = {(N p (s)N q (s') + 
Nq(s')N p (s)))/2 is the steady-state noise correlation ma- 
trix. It turns out that in the regime where the mechan- 
ical oscillators possess very high Q- values, the quantum 
Brownian noise becomes approximately 5-correlated and 
in this limit the noise correlation matrix takes the form 



= Diag[0, 7m (2ii 1 + l)),0, 7m (2n 2 + 1), 

K, K, K, K,]S(s — s'), 



(54) 



with 



3 (ftn/fc B Ti) 



1 



and ri2 = 



jg(ftn/fc B T 2 ) _ j_j \ Thug; Eq ^ q simplifies to 
dsW(s)NW T (s). (55) 



V = 



When the system is stable in the steady state, the covari- 
ance matrix V satisfies a Lyapunov equation [l7j 



zv + VZ J 



-N. 



(56) 
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Once again, when we have the solution for the covari- 
ance matrix, we can compute various non-classical cor- 
relations between the optical and mechanical degrees of 
freedom. In particular, the degree of entanglement be- 
tween different optical and mechanical modes can be eval- 
uated by computing the logarithmic negativity as defined 
in (|2"2"j) . We numerically solve Eq. (f5B|) for the covariance 
matrix V. An example of the numerically calculated log- 
arithmic negativity between various optical and mechan- 
ical modes is presented in Fig. |U 

For evaluating the entanglement between various op- 
tical and mechanical modes we have chosen physical pa- 
rameters accessible in present experiments. Not sur- 
prisingly, the steady-state entanglement is susceptible to 
thermal fluctuations of the environment. A high tem- 
perature of the surrounding reservoirs will result in a 
completely separable state of the optical and mechanical 
modes. One should note that the entanglement generated 
between optical and mechanical modes in the steady state 
is not very large, but, it does not require any quantum 
resources, such as additionally driving the two cavities 
with squeezed light [TH, [l!| . 

Also, it is worth pointing out that with our particular 
choice of parameters we find that an appreciable entan- 
glement appears between various optical and mechanical 
modes only when we operate far away from the regime of 
the red (A = fi) or blue (A = — f2) sideband. Although 
operating in the blue sideband regime is commonly con- 
sidered ideal for generating entanglement between vari- 
ous optical and mechanical modes, the condition that the 
steady state should be stable puts serious restrictions on 
the coupling strength between the mechanical and optical 
modes 5^ . 

A challenging aspect of any scheme involving entan- 
glement generation between macroscopic mechanical sys- 
tems is the actual experimental detection of entangle- 
ment. There are however some recent promising propos- 
als to create and detect quantum correlations in optome- 
chanical settings [ljj|, [22[ • Since it is comparatively easier 
to detect quantum correlations between optical modes, 
as compared to directly detecting quantum entanglement 
between mechanical modes, the essence of these propos- 
als is to swap the nonlocal correlations from the mechan- 
ical modes back to the optical modes. As shown in Fig.Q] 
this can, for instance, be implemented using two auxiliary 
light modes, each initially prepared in classical uncorre- 
cted states. These auxiliary modes can be two modes 
of distant cavities, and the geometry so arranged that 
each entangled mirror couples independently to the two 
modes. The non-local correlations may then be trans- 
ferred from the movable mirrors to the initially uncorre- 
cted auxiliary modes, which may eventually become en- 
tangled. Thus, using standard homodyne measurement 
techniques, the entire correlation matrix of the two op- 



tical auxiliary modes can be reconstructed. A presence 
of non-zero quantum correlations between these optical 
modes will be an indirect signature of non-zero quantum 
correlations between the mechanical modes. 

V. CONCLUSION 

We have discussed in detail the possibility of gen- 
erating nonlocal quantum correlations between optical 
and mechanical modes of two spatially separated cavi- 
ties. Each cavity is assumed to have one fixed and one 
movable mirror and the two cavities arc coupled by an 
optical fiber. Under the Born-Oppenheimcr approxima- 
tion, relying on separating the dynamics into a fast op- 
tical timescale and a slow mechanical timescale, we have 
analytically worked out the dynamics of the two cou- 
pled movable mirrors. Furthermore, within this adiabatic 
regime, we have also presented a full analytical solution of 
the master equation governing the open-system dynam- 
ics of the two movable mirrors. We especially found that 
the interaction mediated via two optical modes entan- 
gles the two distant movable mirrors. Cavity losses were 
taken into account in an effective non-Hermitian model, 
and mirror entanglement was found to be fairly robust 
against such dissipation. Using a complementary model, 
we have studied the two coupled cavities using the quan- 
tum Langevin formalism, by explicitly solving the result- 
ing equations of motion. In the presence of strong driving 
laser fields we have found that the two coupled cavities 
exhibit nonlocal quantum correlations between distant 
optical and mechanical modes. In particular, these op- 
tical and mechanical modes exhibit entanglement in the 
steady state and at finite temperatures. This opens up an 
interesting possibility to study spatially separated mas- 
sive Schrodinger cat states. 

During the completion of this work we became aware 
of a recent proposal discussing the possibility of phonon 
photon entanglement in coupled optomechanical arrays 
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